Ground State Spin Structure of Strongly Interacting Disordered ID Hubbard Model 
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We study the influence of on-site disorder on the magnetic properties of the ground state of the 
infinite U ID Hubbard model. We find that the ground state is not ferromagnetic. This is analyzed 
in terms of the algebraic structure of the spin dependence of the Hamiltonian. A simple explanation 
is derived for the 1/N periodicity in the persistent current for this model. 



The Hubbard model is the simplest model one can 
study to examine the effects of correlations between elec- 
trons in narrow energy bands. The Hamiltonian con- 
sists of a nearest neighbor hopping term and an electron- 
electron repulsion U which acts only when two elec- 
trons are on the same site. The Hubbard model is also 
the canonical model for the study of itinerant ferromag- 
netism. The strong coupling regime is of special impor- 
tance for the study of ferromagnetism, since a theorem by 
Nagaoka jl| states that in the U = oo limit, the GS is fer- 
romagnetic given some connectivity property of the lat- 
tice (which holds in most cases for d > 1). The model is 
solvable in ID and it was shown M that for open bound- 
ary conditions (BC) the ground state (GS) for finite U is 
a singlet, i.e., there can be no ferromagnetism unless one 
postulates explicitly spin- or velocity-dependent forces. 
For infinite U, the GS of all spin sectors are degenerate. 

The problem of the interplay between disorder and in- 
teractions in systems of electrons is challenging and has 
a long history J3| . It is of interest to study the influence 
of disorder on the possibility of forming a ferromagnetic 
GS. In this work we study the spin structure of the GS 
of a disordered U — oo Hubbard model in ID. We find 
that for periodic BC as well as for open BC, for any re- 
alization of on-site disorder, the GS is degenerate, where 
all spin sectors have the same lowest energy, except for 
the fully polarized one which has a higher energy. As a 
by-product of our proof, we find that the GS of an even 
(odd) number of spinlcss fcrmions on a ID ring threaded 
by flux, is minimal when the dimensionlcss flux Q/$o 
equals it (0). This might be of interest for the study of 
persistent currents in disordered interacting ID rings pj . 

Lieb and Mattis g] have considered the one dimen- 
sional clean Hubbard model for any U < oo given hard 
wall (or open) BC. They found that the GS is always a 
singlet (for an even number of spins). When U = oo the 
GS in all the different spin sectors become degenerate. 
Here we present an analysis of the U = oo with periodic 
BC, in the presence of on-site disorder. As will become 
clear later, the different BC change the character of the 
problem, and the periodic BC case give us some insight 
into the higher dimensional variants of the problem. The 
Hamiltonian we consider is thus given by 
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where a irT is the fermionic creation operator on site i with 
spin a, mi = o}+Oif , and the on-site energies Si are drawn 
uniformly between —W/2 and W/2. 

The Hilbert space is composed of a direct product of 
the spatial wave functions, described by a basis composed 
of all the different possibilities of positioning N particles 
on N out of the M different sites. This space is iso- 
morphic to the Hilbert space of TV non interacting spin- 
less fermions on M sites. For finite U, the Hilbert space 
consists of functions with double occupancy as well, and 
these exists only if the two particles on the same site are 
in a singlet configurations. Thus, coupling between the 
spatial functions and spin functions is formed. However, 
in the U — oo case the spin Hilbert space is decomposed 
from the spatial Hilbert space, and its natural basis is 
given by the 2^ orderings of the N spins. 

The action of the hopping terms in the Hamiltonian 
on a wavefunction is changing the spatial distribution of 
the particles, through moving one particle at a time to 
one of its neighboring sites. Since double occupancy is 
restricted, particles can not interchange their order via 
hopping. One gets another invariant of the Hamilto- 
nian, namely, the order of the particles, the importance 
of which will become clear shortly. This holds in the hard 
wall BC case. However, the periodic BC allow particles to 
" bypass" through the boundaries and breaks this symme- 
try. Had the ordering of the spins been conserved, the GS 
would have been independent on the entire spin structure 
of the states. Consequently, in the hard wall BC case, 
all spin configurations GS are degenerate. However, the 
change in the spin ordering due to hopping through the 
boundaries, couples the spin part to the spatial part of 
the wavefunction. Each hopping term in the spatial part 
of the Hamiltonian is replaced by a matrix (in princi- 
ple, 2 N x 2 N ) which characterizes the re-ordering of the 
spins. Most of these matrices, namely, these which are 
not related to the boundary terms, are identity matrices. 
However, some describe the re-ordering due to hopping 
from site 1 to site TV, and vice versa, and are not trivial. 

Let T be the permutation in the spin ordering due to 



hopping of a particle at site 1, to site N. The minimal 
subgroup of the permutations group which contains T 
(and T , which corresponds to hopping in the inverse 
direction) is the cyclic group generated by T. The differ- 
ent spin ordcrings induce a representation of this cyclic 
group. Since the cyclic group is commutative, its irre- 
ducible representations are one dimensional. Thus we 
will find the decomposition of the induced representation 
into irreducible parts, and then diagonalize the spatial 
Hamiltonian for each of the irreducible representations. 
The GS is the lowest of the eigenvalues obtained for each 
representation. Our goal is to determine to which S value 
corresponds the irreducible representation with the low- 
est eigenvalue. 

We note that since T N — /, all the representations 
must be of the form \j(T k ) = exp(2irijk/N). Thus, 
different representations correspond to assigning a phase 
Xj(T) = 2wj/N to the transition through the bound- 
aries. This is equivalent to the adding a flux $ = 2wj/N 
through the ring to the N spinless fermions problem. We 
now prove the following 

Theorem : Let e($) be the GS of N non interacting spin- 
less particles on a (disordered) ring, as a function of the 
flux $. e($) has a minimum at $ = 7r for even number 
of particles (and at $ = for odd N) . 
Proof. The energy as a function of the flux e($), is sym- 
metric in the parameter $ with respect to the points 
<f> = 0, 7r. Therefore, the first derivative e'($) vanishes in 
these points. We now calculate the second derivative at 
these points in order to find which is the maximum and 
which is the minimum. Expansion of e(<£>) in the vicin- 
ity of $ = 7r is given by perturbation theory. The flux 
dependence of the matrix elements is always through the 
expression e**, and therefore one obtains 
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The curvature, i.e., second derivative at $ = -k is thus 
given by 
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At the point $ = n, all the off diagonal elements in the 
Hamiltonian are negative, since the —1 factor of the flux 
which multiplies the boundary hopping matrix elements 
cancels with the (— l)^ -1 factor from the interchange of 
the fermionic operators. Thus the GS eigenvector vq has 
no nodes, i.e., all its elements are of the same sign. Now, 
let us write explicitly a\ and a-i 
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The sum over [ij] is restricted to states i and j which are 
coupled by boundary hopping terms. 

It then follows that the curvature at $ = 7r is positive, 
and therefore it is a local minimum point. Unless there 
is an additional accidental point of vanishing derivative, 
this is also the global minimum. Similar arguments apply 
to the point $ = for odd N. Q.E.D. 

Note that the above theorem holds for off-diagonal dis- 
order as well, as long as all the hopping integrals are of 
the same sign. 

One concludes that for even TV, the irreducible repre- 
sentation which gives the GS is xn/2i an d Xn/2{T) = 
Xn/2(T~ 1 ) = — 1. We now have to find out for which 
S values the induced representation includes Xn/2 in its 
decomposition. There is only one state with S — S z — 
S'max, and thus the induced representation is one dimen- 
sional and irreducible. It is easy to see that it is xo 
(which corresponds to the maximum energy). It will be 
now shown that the representation induced by S = 
always include Xn/2- 

Since the explicit construction of states with definite 
S and S z is non trivial, let us look at the representa- 
tions induced by the sets of states with definite S z only. 
The irreducible representations induced by S z = M, 
which are not induced by S z = M + 1 correspond to 
S = M. According to character theory, the number of 
times an irreducible representation with character x ex- 
ists in the decomposition of a representation (f> is given by 
~W Sfe < fi{k)x*(k : ), where the index k runs through all the 
group elements, and K is the group order. The character 
of a representation is just the trace of the representing 
matrices. The character of the representation induced by 
a set of states is given by the sum over k of the number 
of states in the relevant set which are invariant under the 
fcth element of the group. Accordingly, one obtains that 
the number of occurrences of Xn/2 i n the representation 
induced by S — 1 is given by 



fci 
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N/2 even 
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and the number of occurrences in the representation in- 
duced by S = is bounded by 



fcn < 



(N-l)\ 
(JV/2)!(JV/2)! 



(7) 



It then follows that fco > fci, i.e., there exists at least one 
occurrence of Xn/2 induced by S = states. Therefore, 
the GS is obtained in the S — sector. 

We have shown that for each (even) N, the GS is ob- 
tained at the S = sector, while the lowest energy in the 
Smax sector is higher. 



The above solution is based on the fact that the per- 
mutations induced by the one dimensional boundary hop- 
ping terms spans only the cyclic group, which is commu- 
tative. The representation induced by the spin states is 
highly reducible with respect to this subgroup. In higher 
dimensions, reordering of spins is generated not only by 
boundary terms, and the whole (non commutative) per- 
mutation group is needed to characterize the spin depen- 
dence. It can be shown that the representation induced 
is irreducible with respect to the full non-commutative 
group. Thus, the spin dependence of the Hamiltonian is 
not equivalent to a trivial flux-like correction. 

It is now very easy to consider the effect of magnetic 
flux added to the Hamiltonian. As we have seen, for an 
even number of electrons, the optimal total flux (physical 
flux + flux added by spin configuration) is (in dimension- 
less units) 7r. This indeed is the fictitious flux generated 
by the spins in the GS of the Hamiltonian without an 
external flux. Now, when we turn on the magnetic flux, 
the GS energy will increase. It is therefore favorable for 
the system to produce an inverse fictitious flux to cancel 
out the influence of the magnetic flux. However, as we 
have seen, this fictitious flux comes in quanta of 2ir$$/N , 
where N is the number of particles. Thus magnetic flux 
of integer multiples of 2itQq/N can be completely can- 
celed out by the spin ordering, such that the GS energy 
is exactly as it were in the absence of flux. Therefore, 
as one increases the magnetic flux, the GS energy raises, 
up to the point $ = ir&o/N, where it is favorable for the 
system to produce a negative fictitious flux such that the 
total flux is (in absolute value) less than the magnetic 
flux. We thus have a simple and transparent explanation 
for the flux dependence of the GS energy which is pe- 
riodic with period 2n/N instead of the usual 27r period 



In conclusion, we have shown in this work that the 
GS of spin 1/2 fermions on a ID ring is not polarized in 
the U = oo limit, for any realization of disorder. This 
is accounted for by mapping the spin background of the 
ID problem onto a fictitious flux. In 2D the influence of 
the spin background is non-commutative and therefore 
much more complex. This leads to ferromagnetic GS for 
one hole pi and to disorder induced ferromagnetism for 
higher hole concentration || . It was also shown that the 
energy of spinless fermions is minimized when the ap- 
plied magnetic flux is half the flux period, again, for any 
realization of disorder. This fact is relevant for persistent 
currents calculations. 
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